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We present a new method for the determination of the 
Cabibbo-Kobayashi-Maskawa quark mixing matrix element 
I Vtb I from electroweak loop corrections, in particular those af- 
fecting the process Z — > bb. From a combined analysis of 
results from the LEP, SLC, Tevatron, and neutrino scatter- 
ing experiments we determine |Vt6| — 0.771q24- This is the 
first determination of |Vti,l which is independent of unitarity 
assumptions. 
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I. INTRODUCTION 

The Cabibbo-Kobayashi-Maskawa (CKM) matrix 0] 
describes the relationship between weak and mass eigen- 
states of quarks, assuming that there are three genera- 
tions. By convention, up-type quarks are unmixed such 
that all the mixing is expressible in terms of the 3x3 
unitary matrix V 




Vud Vus Vub 
= I Vcd Vcs Vcb 

Vtd Vts Vtb 




(1) 



where unprimed states denote mass eigenstates and 
primed ones denote weak eigenstates. The unitarity of 
the CKM matrix is an assumption which must be sub- 
jected to experimental verification by independent mea- 
surements of the elements. Non-unitarity of the ma- 
trix would be a clear signature for physics beyond the 
standard model, such as a fourth generation or non- 
universality of the quark couplings. Of particular interest 
in terms of possible new physics is the element Vtb which 
describes the coupling of the two heaviest quarks. 

It is often assumed that |Vtb| « 1 although this ele- 
ment has never been determined without assumptions of 
unitarity. Assuming that there are only three genera- 
tions and unitarity of the CKM matrix yields: 0.9989 < 
\Vtb\ < 0.9993 at the 90% confidence level g. Relax- 
ing the assumption of three generations but maintaining 
that of unitarity yields < |Vffc| < 0.9993 at the 90% 
confidence level Q], while relaxing also that of unitarity 
leaves \Vtb\ unbounded. 

In this paper we describe a new method for the de- 
termination of |Vff,| from electroweak loop corrections, in 
particular to the process Z — > bb. From a combined anal- 
ysis of data from the LEP, SLC, Tevatron, and neutrino 



scattering experiments, we determine for the first time 
the value of |Vtb|, without any assumptions about the 
unitarity of the CKM matrix. The implications for other 
quantities, such as the the top mass, Higgs mass, and the 
strong coupling constant, are discussed. 



II. EFFECTS OF Vtb ON ELECTROWEAK 
RADIATIVE CORRECTIONS 

The precision of the electroweak data from the LEP 
and SLC experiments is sufficient to be sensitive to weak 
loop diagrams involving the top quark. The top quark 
appears in Z vacuum polarisation loops, thereby affecting 
all the Z partial widths, and in the GIM-suppressed ver- 
tex diagrams shown at the one-loop level in Fig. |l| which 
affect the Z partial width, P^^, for the process Z bb. 
From fits to the Z electroweak parameters, the top quark 
mass has been determined to be rrit = 158't\i GeV in 
agreement with the direct measurement from the Teva- 
tron of mt = 175.6 ± 5.5 GeV [| . 

Hitherto the theoretical treatments of weak loop cor- 
rections to Z decay processes have assumed that \Vtb\ — 
1; we relax this assumption. Following the treat- 
ment of Barbieri, Beccaria, Ciafaloni, Curci, and Vicere 
(BBCCV) ||, Tf,b may be written as: 



^bb — 



pRqebRqcdx 1 — 



Ami 
ml 



9bA) 





^^'^ml_ 







(2) 



where is the Fermi constant, mz and mf, are the 
masses of the Z and the 5-quark respectively, and p in- 
cludes the effects of radiative corrections to the Z propa- 
gator. i?QED and i?QCD, which are approximately unity, 
describe the QED and QCD vertex corrections |6|j^ . The 
couplings, gbv and gbA incorporate vertex corrections de- 
scribed by the parameter, r, as follows: 

(3) 
(4) 

where : = ^ | 1 - , / 1 ^1"'" , | . (5) 

We have taken the results of the BBCCV calculations up 
to two-loops and included the effects of Vtb at the one- 
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loop level, such that p is unchanged while r is modified by 
a multiplicative factor of |14hP- In the limit mt ^ rnn, 
T is given by: 



T = -2\VtbVx 



l + -(27-7r^ 



(6) 



where x — G^m^ /87r^V2- In the limit nit <C ran, t is 
given by: 



T^-2\Vtb?x 



We also use the following parameters which are to a 
good approximation experimentally uncorrelated: At = 
—Pt, the average tau polarisation from the 

tau polarisation forward-backward asymmetry ^lr, 
the left-right asymmetry from SLD the QED cou- 
pling constant, a{mz) [||; the strong coupling con- 
stant as{mz) 1^ where the value obtained from the Z 
width is not included; the W masses from LEP II [|, 
CDF §, D0 0, and UA2 flltl, averaged according 
to Ref. [O; the top mass from the Tevatron B; and 



X < 1 + 



144 



(1 



/m^) from z/N scattering measurements by 



311 + 247r2 + 2821ogr + 901ogV CHARM @, CDHS Q, and CCFR (T|. The parame 



-4r (40 + Gtt^ + 15 log r + 18 log^ 



3r2 
100 



(24209 - eOOOTT^ - 45420 log r 



-18000 log r) 



(7) 



where r — mf/mj^. In the intermediate region {mt ~ 
mn) T is described by a polynomial parametrisation of 
the full BBCCV calculation, as a function of mt/rriH, 
multiplied by a factor of |Vtfcp. At the two- loop level 
some diagrams are of the order \Vtb\'^- Our treatment of 
these to order |VtbP is justified by their relatively small ble |IV 
contribution and the current sensitivity of the experimen- 
tal data, as will be seen below. 



III. DETERMINATION OF |Vtb| FROM A FIT TO 
ELECTROWEAK DATA 



ter values and their errors are shown in Table [II . 

For given values of mz, nit, thh, cts, Qf, and |Vth| 
our modified version of ZFITTER provides predictions 
for all of the parameters shown in Tables |, ||, and III. 
These predictions and the corresponding measured quan- 
tities, together with their associated errors and correla- 
tion coefficients, are used to construct a chisquare prob- 
ability ^^(w^, mt, m^, Q!s, Q;, Vtb). The minimum of the 
chisquare is then determined numerically. As a technical 
cross-check, we use the same input parameters as those 
of Ward set |Vtf,| = 1, and successfully reproduce the 
results for mz, mt, mn, ctg, and a. 

The results of the fit with \Vtb\ free are shown in Ta- 
We determine \Vtb\ = 0.77to;24. This value is con- 
sistent with the unitarity prediction of |Vtb| « 0.9991 to 
within approximately one standard deviation. For com- 
parison. Table pv| also includes the results with |Vtb| fixed 
to 0.9991. In both fits, the probability is consistent 
with expectations given the number of degrees of free- 
dom. 



The BBCCV corrections are incorporated in the ZFIT- 
TER program which is used by the LEP/SLC elec- 
troweak working groups to derive parameters such as mt 
and mn from Z data ||^ . We make modest modifications 
to ZFITTER to allow for the effects of Vtb described in 
the previous section. 

The Z parameters, from a combined fit to LEP/SLC 
data, which we use as input are [^: the Z mass, mz', 
the Z width, Tz', the hadronic pole cross-section, crJJ; 
Ri = rhad/r« where Fhad is the hadronic partial Z 
width and Tu is the leptonic partial width, assuming lep- 
ton universality; and the leptonic pole forward-backward 
charge asymmetry assuming lepton universality, Apg. 
The parameter values, their errors, and their correlation 
coefficients are shown in Table |. 

The parameters pertaining to b and c quarks which 
we use are §: i?^ = T.-jTh^d] i?° = Tee/rhad; ^fb 
and Ap^, the forward-backward charge asymmetries at 
the Z pole for b and c quarks respectively; and Af (f — 
b,c) where Af = 2gfvgfA/ {gj^y + <7/a)- parameter 
values, their errors, and their correlation coefficients are 
shown in Table ||. 



IV. DISCUSSION 

The fitted values of mz and a{mz)~^ are insensitive 
to I Vtb I as expected, as shown by the correlation coeffi- 
cients from the fit with |Vf6| free in Table 0. The anti- 
correlations of |Vtf,| with mt and m^, shown in Fig. |^(a) 
and |(b) and in Table 0, have only a weak effect on the 
determinations of mt and mn since the Tevatron mea- 
surement of mt and the vacuum polarisation contribution 
to the Z width constrain mt and mn independent of |Vtb|. 
Nonetheless, allowing |Vth| to float increases the fitted 
values of mt by 1.5 GeV and of mn by approximately 
30GeV. Allowing |Vtb| to float decreases the fitted value 
of as by approximately 0.7 standard deviations, due to 
the fairly strong correlation of these two quantities, as 
shown in Fig. |^(c) and in Tabic 0. 

To assess the future sensitivity of this technique for 
determining |Vth| we reduce the error by a factor of two 
on each of the input parameters in turn, without chang- 
ing the errors on the other parameters. The only pa- 
rameters which cause A|Vtb|/|Vtb| to change by a relative 
amount of more than 10%, from the original value of 
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A|Vtfc|/|yt6| « 26.5%, are T^, 7?°, and a,(mz), which 
yield uncertainties of lS.\Vtb\/\Vtb\ « 22.5%, 19.9%, and 
22.4% respectively. Factors of somewhat less than two 
may be expected from the final analyses of the LEP I 
and SLC data samples. We estimate that the error on 
|Vtb| from the final Z samples of the LEP and SLC ex- 
periments will be approximately 20%. 

Recently CDF presented preliminary results of an anal- 
ysis of their tt event samples. However, they necessarily 
assumed that there are only three generations and that 
|Vtdp + |VfsP + I Vtbp = 1, as required by unitarity, to ex- 
tract iFtbl^son ^ o.99±0.15 ge). Ultimately, the measure- 
ment of the single top quark production rate at hadron 
colliders should be sensitive to \Vtb\ without requiring 
such assumptions of unitarity. The estimated sensitivity 
at the end of the Tevatron Run II for (5|Vifc|/|yt6| is 12%- 
19%, depending on the uncertainty of the gluon structure 
functions [f?!. 



V. SUMMARY 

We describe a new technique for the determination of 
the CKM matrix element |Vtb| using loop corrections to 
electroweak processes, without using any unitarity con- 
straints. From a combined analysis of data from the LEP, 
SLC, Tevatron and neutrino scattering experiments we 
determine |Vtb| = 0.77lg'24 where the error includes the 
experimental errors and uncertainties on the top mass, 
the Higgs mass, and the strong coupling constant. The 
wider implications of this measurement are discussed 
elsewhere [ITHl. 
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FIG. 2. Variation of |Vt6| with (a) mt, (b) logio(mH/GeV), 
and (c) a,. The point with error bars denotes the result of the 
fit allowing for all errors and correlations. The solid line shows 
the dependence of |Vt(,| on the ordinate variable; the dashed 
lines correspond to the 68% confidence level. The hatched line 
of (b) shows the low vhh region excluded by direct searches 
at LEP II. 



TABLE I 



Measured values and the correlation matrix for 
Ri, and J^^L from a combined fit to LEP and 



SLC data. 



Measured 



Correlation coefficient 





value 


mz 


Tz 


<rl 


Re 


^FB 


mz (GeV) 91.1867 ± 0.0020 


1.00 


0.05 


-0.01 


-0.02 


0.06 


Vz (GeV) 


2.4948 ± 0.0025 


0.05 


1.00 


-0.16 


0.00 


0.00 


al (nb) 


41.486 ±0.053 


-0.01 


-0.16 


1.00 


0.14 


0.00 


Re 


20.775 ±0.027 


-0.02 


0.00 


0.14 


1.00 


0.01 


^FB 


0.0171 ± 0.0010 


0.06 


0.00 


0.00 


0.01 


1.00 



TABLE II. Measured values and the correlation matrix for 

from a combined fit to LEP 



O.b 

and SLC data 



R'L i?" A 



Ap'g, Ab, and Ac 



Measured 



Correlation coefficient 



Rb 
R^ 

AO,b 

^FB 
.0,c 
^FB 

At 
Ac 



value 


Rb 


Rc 


^FB 


^FB 


Ab 


Ac 


0.2170 ± 0.0009 


1.00 


-0.20 


-0.03 


0.01 


-0.03 


0.02 


0.1734 ± 0.0048 


-0.20 


1.00 


0.03 


-0.07 


0.04 


-0.04 


0.0984 ± 0.0024 


-0.03 


0.03 


1.00 


0.13 


0.03 


0.02 


0.0741 ± 0.0048 


0.01 


-0.07 


0.13 


1.00 


0.00 


0.07 


0.900 ±0.050 


-0.03 


0.04 


0.03 


0.00 


1.00 


0.08 


0.650 ±0.058 


0.02 


-0.04 


0.02 


0.07 


0.08 


1.00 



TABLE III. MeEisured values of uncorrelated parameters 
used in our fits. 

Measured value 



Ar 
Ac 

a^^{rnz) 
as{mz) 
mw (GeV) 
mt (GeV) 

1 - n}\ylivi^ 



0.1410 ±0.0064 
0.1399 ±0.0073 

0.1547 ±0.0032 
128.896 ± 0.090 
0.118 ±0.003 
80.400 ± 0.075 

175.6 ±5.5 
0.2251 ± 0.0037 



TABLE IV. Results of the fit for mz, mt, log^(,(mi//GeV), 
as(mz), a{mz)''^, and \Vtb\ (second colurrm). For compar- 
ison, the third column shows the results of the fit with the 
constraint of |Vti,| = 1. P denotes the probability of obtain- 
ing a reduced chisquare greater than that from the fit. 





\Vtb\ free 


\Vtb\ fixed 


mz (GeV) 


91.1866 ± 0.0020 


91.1866 ± 0.0020 


mt (GeV) 


174.2 ± 5.4 


172.7 ±5.2 


logio('TiH/GeV) 


9 1 cr+O.SO 
^•J-J_0.39 


2.04l°j? 


as(mz) 


0.1171 ± 0.0025 


0.1188 ±0.0021 


a~^(mz) 


128.913 ± 0.092 


128.905 ±0.091 


\Vtb\ 


^- ' ' -0.24 


0.9991 (fixed) 


xl = xVd-o.f 


15.1/(19-6) 


16.7/(19 - 5) 


> 5S) (%) 


30 


27 



TABLE V. Correlation eocffieieuts from the fit for mz. m,t, 



logi()("'-r//Ge\'), a„ 




a(itiz 


) \. ami 










mz 


mt 


logio 


imn) 


Ob 




1^*6 1 


mz 


1.00 


0.01 




0.04 - 


-0.01 


0.01 


0.01 


mt 


0.01 


1.00 




0.65 - 


-0.03 


0.16 


-0.15 


logio(n^ii) 


0.04 


0.65 




1.00 


0.04 


0.65 


-0.23 


as 


-0.01 


-0.03 




0.04 


1.00 


0.03 


0.52 




0.01 


0.16 




0.65 


0.03 


1.00 


-0.08 




0.01 


-0.15 




-0.23 


0.52 


-0.08 


1.00 
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